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Abstract 


The hybrid solid oxide fuel cell (SOFC) and gas turbine (GT) system is a promising concept in the future power generation for its high- 
performance and low-emission. The dynamic model for the hybrid system of integrated SOFC and recuperative GT with air reheating component 
is presented in this work. A dynamic model was put forward based on the conservation equations of mass, energy and force through the whole 
plant, with specific source terms in different types of components. The SOFC was modeled on the basis of the Exponential Decay function and the 
Exponential Associate function, which describe the characteristics of the parameters distribution within the SOFC. A cubic curve was employed 
to denote the compressor pressure characteristics. In the turbine model, the relation between the work done and the inlet condition of turbine was 
determined according to the turbine nozzle work characteristics. The developed system model was programmed and implemented in the simulation 
tool Aspen Custom Modeler. The current density of SOFC was selected as disturbance variable during the dynamic simulation using the developed 
dynamic model. The responses of the SOFC air inlet temperature, SOFC outlet temperature, and turbine inlet temperature, the output voltage, and 
the gas species molar fractions at the outlet of SOFC were presented. The obtained results show that the presented dynamic model can be able to 


simulate the system dynamic track reasonably. 
© 2006 Elsevier B.V. All rights reserved. 


Keywords: Solid oxide fuel cell; Hybrid system; Dynamic model 


1. Introduction 


The hybrid solid oxide fuel cell (SOFC) and gas turbine (GT) 
system is a promising concept in the future power generation 
for its high-performance and low-emission. The interest in the 
hybrid system has been increasing with the improvement of GT 
technologies and maturity of modular fuel cell in recent years, 
leading to the fact that the high-performance hybrid cycle can 
be applied in small-capacity power systems. A research project 
Vision 21 program [1] was issued by the United States Depart- 
ment of Energy in 1997 for conceptual feasibility studies of 
high-efficiency fossil power plants, which aimed to develop the 
core modules for a fleet of fuel-flexible, multi-product energy 
plants with the power efficiency higher than 75% for the gas- 
fueled. The integration of fuel cell and GT hybrid systems was 
considered to be the key technology to achieve the Vision 21st 
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goals due to its high-performance and low-emission. Siemens- 
Westinghouse Power Corporation has issued several conceptual 
design of pressurized SOFC/GT hybrid system [2]. The first 
demonstration project of a pressurized SOFC/GT hybrid sys- 
tem provided by Siemens- Westinghouse Power Corporation was 
carried out successfully at the University of California-Irvine. 
Furthermore, Allison Engine Company evaluated a pressurized 
SOFC coupled with conventional GT technology plant. The sys- 
tem was designed to operate at a pressure of 7 atm and the 
process showed that 67% efficiency was reached for power 
generation. An efficiency of 70% was deemed achievable with 
the improvement in component design. The cost of electricity 
(COE) was predicted to be comparable to the alternatives of 
current. Besides, NOx emission was less than 1 ppm [3]. In addi- 
tion, McDermott Technology developed a conceptual design of 
a high efficiency power plant system which joined the planar 
SOFC technology with microturbine technology together in a 
combined cycle. The system was operated at atmospheric condi- 
tions while the output of the combined cycle for the power plant 
achieves 700 kW. The results indicated that 70% efficiency was 
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Nomenclature 


ao, 41, a2 the fitting function coefficient of specific heat 


capacity (Jmol~! K7!) 

area (m2) 

the parameter of exponential fitting function 
coefficient of discharge 

specific heat capacity at constant pressure 
(Ikg~! K) 

activation energy (J mol~!) 

theoretical Nernst potential (V) 

molar ratio of methane to water 

the reforming reaction balance constant 
Faraday constant (=96,485 C mol~!), force (N) 
specific enthalpy (Jkg~!) 

change of specific enthalpy (Jkg~!) 
parameter of formula (9) 

current density (A m~?) 

diffusion limiting current density (A m~?) 
exchange current density (A m~?) 

specific kinetic energy (Jkg~!) 

reforming chemical reaction constant 
(=4274 mol s7! m~? bar!) 

shift gas—water chemical reaction constant 
(=1.2 x 104 mol m~? s7!) 

equilibrium constant 

mass flow (kg s7!) 

molecular mass (kg mol`!) 

molar weight (mol) 

molar flow (mol s~!) 

electrons transferred per reaction 
operating pressure (Pa) 


Preference reference pressure (Pa) 


standard pressure (Pa) 
chemical reaction heat (J mol~!) 


Qexchange exchanged heat (J s7!) 


generated heat (J s7!) 

molar universal gas constant 

(=8.314J mol~! K7!), chemical reaction rate 
(mol s~!) 

electrolyte ohmic resistance ({2) 

mass universal gas constant (Jkg~! K7!) 
mass source (kg s7!) 

species molar source (mol s7!) 

heat source (J s7!) 

work done (W) 

time (s) 

fitting function exponential parameter 
temperature (K) 

impeller tip velocity (m s7!) 

voltage (V) 

fluid velocity (m s7!) 

volume (m?) 

parameter of formula (9) 

gaseous molar fraction 


y fitting function 
Zz dimensionless spatial coordinate of gas flow 
direction (z € [0,1]) 


Greek letters 

ratio of specific heat 

thickness of housing for control volume (m) 
heat exchanger effectiveness 

turbine flow factor 

stoichiometric reaction coefficient 

splitter diffluent factor 

turbine work done factor 

density (kg m7?) 

mass flow coefficient 

X corrected coefficient of Eq. (3) 

w, Wo isentropic work coefficient 


SDd IlIno or 


Subscripts 
active chemical active 
av average 


cold cold side control volume 
CH4 methane 


CO carbon monoxide 

CO2 carbon dioxide 

el electrochemical reaction 
hot hot side control volume 


H2 hydrogen 
H20 water (gas) 


i index of the control volume 
imp impeller 

isen isentropic 

j index of gaseous species 

k index of anabranch flow 
nozzle turbine nozzle 

N2 nitrogen 

O2 oxygen 

ref reforming reaction 

rot in rotating equipment 


shift shift reaction 
1-6 index of the chemical reactions (R1)—(R6) 


Superscripts 

add addition 

anode anode side control volume 
cathode cathode side control volume 


in inlet 
out outlet 
s solid 


possible and the COE was comparable to present day alternatives 
[4]. 

The GT is driven by employing the rejected thermal energy 
and residual fuel from a fuel cell in the SOFC/GT hybrid system. 
Fig. 1 shows a cycle scheme of SOFC/GT hybrid system with 
recuperative heat exchanger (RHE) as air reheater. Some models 
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Fig. 1. Scheme of a SOFC/GT hybrid system with RHE air reheater cycle. 


adapted to the cycle of hybrid system were presented in Refs. 
[5-8]. Most of them were developed on the basis of a steady 
point or a specific product for the balance of plant (BoP) of 
the hybrid system. The components of compressor and turbine 
dynamic work performance characteristics were not considered 
enough in the system model, as well as the spatial variation 
effect within the SOFC, which was not taken into account in the 
lumped SOFC model [9]. 

The aim of present work focuses on the development of a 
dynamic model for the SOFC/GT hybrid system. In Section 2 of 
this paper, we build the models based on the mass balance differ- 
ential equation, energy balance differential equation and force 
balance equation for each of the components. Section 3 presents 
the dynamic simulation carried out using the developed model, 
which is programmed and implemented in the simulation tool 
Aspen Custom Modeler [10] allowing users to develop custom- 
made model. Then the responses of some parameters such as 
temperature, fuel cell output voltage and species molar fraction 
are analyzed while the current density of the SOFC was assumed 
as the disturbance variable during the simulation. A summary 
of the work is given in Section 4. 


2. Mathematical model 


Firstly, model of each component is built for the hybrid sys- 
tem. The most important assumptions used in present study are 
listed as follows: (1) one-dimensional treatment of flowpaths. 
All fluid properties are assumed to be uniform over any cross 
section because the rate of fluid properties variations perpendic- 
ular to the streamline is small compared with the rate of changes 
along the streamline; (2) the friction loss along the streamline 
is neglected; (3) both air and fuel are assumed as the ideal gas, 
which makes them satisfy the ideal gas state equation; (4) the dif- 
ference of the fluid momentum between the inlet and the outlet 
for each component is neglected. 

In order to make a process model be easily adapted to varia- 
tion of the process, it is necessary to build a generic model. In 


this way, system model should be easily extended or modified 
such as layout schemes, components size. The control volume 
is defined for the place, where the gas state is varied within the 
component along the flowpaths. A general model for each of 
control volume is to be proposed, and then the specific source 
terms are to be determined. 

The mass balance, energy balance and gas species molar 
balance should be satisfied within the control volume. The cor- 
responding conservation equations are written as: 

For the mass balance: 


doc 
yi = 
dt 


where Sm is the term of mass source for the control volume. 
For the energy balance: 


„in + out 
Mi — Mi + Smi 


(1) 


Lo Vic’, + V;,0e% (Cott — Rov Tha 
PLN pee YiPi Mpi — Agi H- a 


= hP aP + K) — ghp" + KM) + Sgi + Swi 2) 
where Sq is the heat source and Sw is the work done for the 
control volume. The gas mass constant Rg is obtained accord- 
ing to the composition of the gas mixture. The specific heat at 
constant pressure is assumed quadratic relation with tempera- 
ture. The term pf V? ji represent the housing thermal capacity 
of the control volume. In this paper, the material of the control 
volume housing is assumed to be Ferro and the control volume 
is assumed to be tetragonal volume. Then the solid volume of 


housing can be estimated according to the follow equation: 


(3) 


where 6 denotes the thickness of the control volume housing. 
The coefficient x is introduced for considering the control vol- 
ume that may be includes solid component. The fuel cell stack 
thermal capacity is calculated according to the single-unit cell 
geometry parameters, the total active area, and the heat proper- 
ties of electrode, electrolyte and connector. 


Vs = 6x v?” 
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For each of gas species molar balance equations 


dn?” 
IJ ein in „out ont ig (4) 
a ~ ni Xij ni Xij nij 


where Sn is species molar source for the control volume. Molar 
weight nj satisfies ideal gas equation of state. The force balance 
is assumed for each of the control volume, then, we get 

pi Ay = preva” + Frrictioni + Froti + Fareachange i =0 (5) 
The friction force is negligible according to the assumption (2). 
Fareachange i 18 the force produced by the pressure effect in case 
of the variation of the cross section area along the flow direc- 
tion within the control volume. The flowing area is assumed as 
constant within the control volume in this study. Then Eq. (5) is 
simplified as follows 


Aip? — 


where Frot is the force on the fluid in the control volume. In this 
paper, it was calculated only in compressor and turbine compo- 
nents in the hybrid system. 


Pe") + Froti = 0 (6) 


2.1. Compressor model 


In order to simulate the compressor behavior, relations have 
to be found among the source terms (forces on, heat transport to, 
and work done on the fluid), inlet condition, mass flow, speed 
and geometry of compressor. The fluid isentropic enthalpy rise 
is first determined according to the flow coefficient. The flow 
coefficient is defined as: 

yin 


oe 7) 


u 


The fluid isentropic enthalpy rise is non-dimensionalized 
through division by the square impeller tip velocity: 


=; (8) 


A cubic relation between the fluid isentropic enthalpy rise and 
the flow coefficient is adopted according to the work of Moore 
[11]. Fig. 2 shows one such curve used in this paper. The formula 
of the cubic curve is 


7 3(¢ 1/¢ i 
vø=w+nji+ ($ 1) (4 1) (9) 


where the parameters Wo, H, and W are illustrated in Fig. 2. 
The isentropic enthalpy rise Ahisen is obtained through Eqs. (8) 
and (9). Then the force Frot on the gas in the impeller for force 
balance Eq. (6) are calculated from Ahisen using 


i i H(y-V) 
p™ + Fyot/Aimp av — T™ + Ahigsen/ Cp es 10 
p” m Tin (10) 
Then 
Sw = rh?” Ahisen (11) 
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Fig. 2. Typical compressor enthalpy rises [11]. 


2.2. Turbine model 


Like the approach taken to formulate and construct the com- 
pressor model, turbine force and work are also determined 
from a set of turbine performance characteristics. To describe 
the turbine performance characteristics, two parameters are 
defined—turbine flow factor, 0, and turbine work done factor, 
T, 


my T 


9 = = (12) 
—Ahisen 


A typical turbine performance curve [12] is drawn in Fig. 3. In 
this paper, the turbine nozzle area is set to be a variable, and then 
the expansion ratio of the turbine is modified by changing the 
nozzle area. Assuming the expansion ratio at the nozzle choke 
condition is known, we can obtain the value of choke based on 
the assumption of isentropic process for the turbine. The choke 
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Fig. 3. Typical turbine performance curves [12]. 
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flow of the turbine nozzle is calculated by the following equation: 


m = CJA in| Y 2 
choke dAnozzle P RTin y+1 


(y+ D/(y-1) 
) (14) 


The turbine flow factor is determined from the inlet condition 
according to Eq. (12). Then, the work done factor z is calculated 
based on the turbine performance curve. The source of output 
power becomes: 


Sw = am™y™ (15) 


The force produced by gas fluid on the rotor is expressed with 
Eq. (10), in which the enthalpy change term is calculated from 
the work done factor. 


2.3. SOFC model 


A planar SOFC developed by International Energy Agency 
(IEA) [13,14] was used to be combined with GT in this study. 
The electrochemical reaction formula for the planar SOFC is 
(1/2)02 + H2 > H20 (R1) 

For the internal reforming of presently used SOFC, the 
reforming reaction and gas-shift reaction are included in the 
model: 


Steam-reforming : CH4 +H20 — CO + 3H? (R2) 


Shifting : CO + H20 = CO2 + H2 (R3) 

The geometry parameters and the design operating condi- 
tion of the SOFC were given in Ref. [14]. A mathematic model 
was presented in [9], which employed two fitting functions—the 
Exponential Decay function and the Exponential Associate func- 
tion to describe the characteristics of the gaseous molar fraction 
and temperature distribution in the flow direction in the SOFC. 
The fitting functions were expressed as: 


The reforming reaction rate is 


1 
Ex 
Rref = Aretkreform reform P f XCH4 EXP (- a) dz (18) 
0 


where xcy,, T are the functions of z satisfying the formulations 
(16). The water—gas shift reaction is given as follows 


1 
Rsnift = AshiftKshift P 1 (xcoxmo — KsniftXco.XH,) dz (19) 
0 


where Kgnift is obtained using the method in Ref. [15]. The 
sources of the mass, the species molar, and the energy balance 
equations are calculated, respectively with 


gue = —Mo,A0,,1 Rel (20) 
geathode = Mo,A07,1 Rel (21) 
Sande — A j1 Rel + AjrRret + Aj,3 Rshifi 22) 
el = Aji Rel (23) 
anode th anode 

anode __ p 

Sq ~~ Canode;j, anode 4 (cathode cathode Ogen (24) 
P p 

, ce athode a cathode 

See = a Osen (25) 


Canode m anode + Ccathode m cathode 
P P 


where j denotes H2, H20, CO, CO2, CH4 for anode side and O2 
for cathode side. And the released heat of chemical reaction is 
calculated according to the following formula: 


Q gen = (qel Rei — Ui Aactive) + Gref Rref + qshift Rshift (26) 


The Nernst potential is calculated as a function of current 
density, operating pressure, temperature and the gaseous molar 
fraction. The irreversibility in the voltage drop is contributed by 
three kinds of polarization losses, i.e. activation, concentration 
and ohmic overpotentials. The Nernst equation for the SOFC is: 


yo + (Ape &/) + A2e7€/2) for the Exponential Decay function 


— yo + [Ao — e7 &/) + A2(1 — e~ &/2))] for the Exponential Associate function 


where y denotes the gaseous molar fraction and the temperature, 
which is a function of the dimensionless spatial coordinate of 
gas flow direction z. The parameters of ti, t2, yo, A1, A2 are 
determined according to the method from Ref. [9]. 

In order to model the sources of the energy balance equation, 
mass balance equation, and gas species balance equations, the 
rate of the electrochemical reaction, the reforming reaction and 
the water—gas shift reaction should be known. The electrochem- 
ical reaction rate is defined as 

i 
Rey = Aactive > (17) 
neF 


2 
RT RT x 0 
i= In K(T) In moe (27) 
neF 2ne F Xf X02 
(16) 
where the assumption p2"¢? = p°thode = p is proposed. The out- 


put voltage can be obtained through Eq. (27) subtracting the 
electrochemical losses 


1 
ü=] 
0 


. 2RT 
En —iR.(T) — 
n 


where the second term on the right hand side of the above equa- 
tion represents ohmic loss, the third term is the activation loss 
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at the anode side, the fourth term is the activation loss at the 
cathode side, and the last term is the concentration loss. The 
exchange current density is modeled in the following [16]: 


anode anode anode 
janode _ pinode P XH) Pp XH,O exp EN 
0 = = 
Preference Preference RT 


(29) 
0.25 
hod hod 
jonthode = oathode pono’ x XO exp | — EX cia (30) 
Preference RT 


2.4. Pre-reformer model 


The reforming reaction and water—gas shift reaction are 
included in the pre-reformer. The formulas of the reforming 
reaction rate and the shift reaction rate are obtained according to 
Eqs. (18) and (19) with no integration. The source term of mass 
balance equation equals to the needed water, which is added to 
ensure the reforming reaction proceeding. Then 


„inin 
n XCH4 _ in in 
m = —— —-n XO Mp,0 (31) 


where fgw is the molar ratio of methane to water. For the species 
molar fraction balance equations, there have 


Sn, j = Aj2Rret + À j,3 Rshift (32) 


where j denotes H2, CO, CO2, CH4. The added water should be 
considered in the water conservation equation: 


in in 
n XCH H4 „inin 
fe = XH,O 
W 
(33) 


Sn, H20 = ÀH20,2 Rref + AH0,3 Rshift + 


The additional heat is needed because of the endothermic 
reforming reaction. Therefore, 


Sq = dref Rref + qshift Rshitt + Smhi (34) 
2.5. Combustor model 


In the combustor the combustion of the exhaust hydrogen and 
carbon monoxide from fuel cell should be taken into account. 
The additional fuel is added to control the turbine inlet temper- 
ature. Three oxidation reactions exist in the combustor: 


CH4 + 202 = 2H20 + CO2 (R4) 
CO + 0.502 = CO2 (R5) 
H2 + 0.502 = H20 (R6) 


The source of mass balance equation is exactly equal to the fuel 
added. Here we assume the combustor efficiency to be constant 
and equal to 100%. Then we obtain: 


— 3., gin in fuel „fuel „pin „in 
Snj = Misra NCH, +H, xCH,) + Ais@ico*co 


+ AES SES) + Aoi af + ARE 65) 


The heat added to the energy equation is given as 


„in in fuel „fuel - fuel „fuel 
Sq = qCH4 (CH, X CH, HÊCH4XCH,) + qcolneoxto +1¢0 co) 


+ guy (iy, xy + AS HS) (36) 


2.6. Heat exchanger model 


The model of heat exchanger consists of two control 
volumes—a cold flow and a hot flow. The heat exchange between 
the two flows is calculated using a given recuperator effective- 
ness [17]: 


in 
Qexchange = emin(CHroiaeoias C phot! hor)(Ty hot — Koa) (37) 


In this work e is set to be 0.85. The source of energy equation 
then becomes 


Sq cold = Qexchange (38) 


Sq hot = — Qexchange (39) 


2.7. Mixer model 


For the case of two streams or multi-stream confluence, the 
maximal mass stream is defined as mainstream, which satis- 
fies the balance equation. The anabranch flows are specified as 
the sources of mainstream. The force on the mainstream by the 
anabranch flow is assumed to be zero in the mainstream flow 
direction. Then 


Sm = X rn (40) 
k 
= 2 me xe) (41) 


Sq = DD iba he; (42) 


where the subscripts k, i denote the index of anabranch and gas 
species, respectively. 


2.8. Splitter model 


Like the confluence mixture model, a maximal mass 
anabranch flow is defined as the mainstream. A variable named 
diffluence factor is first proposed to control the corresponding 
anabranch mass flow, which is defined as 


mou 


& = 


(43) 


min 


Then the sources for balance equations are determined according 
to the following formulas: 


Sm = So (~ern) (44) 
k 
Sai = 25 thg xE) (45) 
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Fig. 4. Disturbance course of current density. 
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where the subscripts k, i denote the index of anabranch and gas 
species, respectively. 
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3. Dynamic simulation 


A planar SOFC presented by IEA is adopted as the simulation 
case in this work. The geometry parameters of the single-unit 
cell are given as follows: active area 5.42 mm x 100 mm, anode 
thickness 50 um, cathode thickness 50 um, electrolyte thickness 
150 wm, channel height 1 mm, rib width 2.42 mm, and bipolar 
plate thickness 2.5mm. The properties of electric conductiv- 
ity for the electrodes and the electrolyte were specified in Ref. 
[13]. The mass density and heat capacity of anode, cathode, 
electrolyte, and connector was assumed to be 6600 kg m~? and 
400 J (kg K)~!, respectively according to IEA report [13]. The 
fuel cell stage is composed of single-unit cells. The designed 
parameters of the proposed system are given in Table 1. 

The simulation is carried out using the proposed model in 
the simulation tool Aspen Custom Modeler, which allows users 
to build custom-made mathematical model. All the component 
models are programmed with the modeling language provided 
by Aspen Custom Modeler programmer environment. 

The current density of SOFC is selected as the disturbance 
during the dynamic simulation. Fig. 4 shows a specific distur- 
bance characteristic of the current density. The fuel mass flow 
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Fig. 5. Response of temperature and output voltage: (a) SOFC outlet temperature; (b) SOFC air inlet temperature; (c) turbine inlet temperature; (d) SOFC output 
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Fig. 6. Response of the gases species molar fraction after the fuel cell: (a) H2; (b) CO2; (c) CO; (d) H20; (e) 02; (f) No. 


and air mass flow of fuel cell are assumed to be fixed during the 
simulation. Steady simulation is first implemented at the design 
point. Then the dynamic mode is set and the simulation is run 
at the design operating condition, in which the current density 
equals to 3000 A m~?. The current density jumps up to the level 
of 3100 A m~? firstly and then down to 2900 A m~ after a few 
hours according to Fig. 4. 


In Fig. 5 the responses of temperature are presented, including 
the temperature of SOFC air inlet (a), SOFC outlet (b) and the 
turbine inlet (c). The polarization losses of SOFC are positively 
correlated with current density [18]. Therefore, the released heat 
increases as current density rises. Then as shown in Fig. 5(b), the 
response curve of SOFC outlet temperature changes as the same 
direction with current density disturbance. However, the fuel 
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Table 1 
Assumed system components performance 
Compressor Isentropic efficiency 85% 
Inlet air pressure (Pa) 101,325 
Inlet air temperature (K) 298 
Inlet air composition 02: 21%, N2: 79% 
Shaft speed (rpm) 46,000 
Pressure ratio 5.157 
Turbine Isentropic efficiency 85% 
Inlet gas temperature (K) 1200 
Expansion ratio 4.55 
SOFC Fuel utilization 85% 
Air utilization 25% 
Current density (A m~7) 3000 
Inlet temperature (K) 1173 
Combustor Efficiency 100% 
Heat exchanger Effectiveness 85% 
Pre-reformer Steam/methane ratio 2.5 
Operating temperature (K) 1173 
Inlet fuel composition 100% CH4 


utilization decreases with increasing of current density when 
the fuel mass flow is set to be fixed. In the other words, the 
residual fuel from the exhaust of SOFC that is combusted in 
the combustor to raise the mixture gas temperature decreases as 
current density rise. So the turbine inlet temperature and the air 
inlet temperature of SOFC are negatively correlated with current 
density. It can be seen that the track of above three temperature 
variations due to the current density disturbance is simulated 
reasonably. Fig. 5(d) shows the response of the output voltage 
of the SOFC, which increases with the decrease of the current 
density. 

Fig. 6 shows the responses of the SOFC outlet gaseous molar 
fraction. The responses of each of the species molar fraction are 
much faster than that of temperature. The steady value of the 
carbon dioxide, water and nitrogen molar fractions are in pro- 
portion to the current density while the variations of hydrogen, 
oxygen and carbon monoxide molar fraction go in the opposite 
change direction of the current density. From Fig. 6(b) and (c), 
we can find that the fluctuating direction of the carbon diox- 
ide molar fraction around the steady value at the jump point is 
inversed with that of carbon dioxide molar fraction. This is due 
to the fact that the shift reaction is closely related with temper- 
ature. The shift reaction shows more negative direction as the 
operating temperature of SOFC increases. So the molar frac- 
tions of reactant species CO increase and the product species 
CO, decrease if the temperature rises. 


4. Conclusion 


A dynamic model of the SOFC combined with recuperative 
gas turbine hybrid system has been developed. Dynamic simu- 


lation for the proposed hybrid system was implemented based 
on the developed model, which was programmed via the simu- 
lation tool Aspen Custom Modeler. The current density of the 
SOFC was selected as disturbance variable in the simulation. 
The fuel and air mass flows were fixed and the fuel utilization 
and air utilization were specified as free variables during the sim- 
ulation. Dynamic simulation results show that the SOFC outlet 
temperature was positively related with current density while 
was inversed for SOFC air inlet temperature and turbine inlet 
temperature. The response time constant of the species molar 
fraction is much smaller than that of temperature. Simulation 
results show that the developed model can follow the track of 
the disturbance very well. 
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